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Abstract. We present new results on the kinematics of the jet HH 110. New proper motion measurements have 
been calculated from [SII] CCD images obtained with a time baseline of nearly fifteen years. HH 110 proper 
motions show a strong asymmetry with respect to the outflow axis, with a general trend of pointing towards 
the west of the axis direction. Spatial velocities have been obtained by combining the proper motions and radial 
velocities from Fabry-Perot data. Velocities decrease by a factor ~ 3 over a distance of ~ 10 18 cm, much shorter 
than the distances expected for the braking caused by the jet /environment interaction. Our results show evidence 
of an anomalously strong interaction between the outflow and the surrounding environment, and are compatible 
with the scenario in which HH 110 emerges from a deflection in a jet/cloud collision. 
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1. Introduction In contrast to the optical images, the HH 110 jet 

is nearly straight in near infrared H2 2.12 am. images. 
The HH 110 jet (which is found in the Orion B clou d com- Furthermore, the H 2 emission appears shifted westward 
plex) was discovered by Reipurth & Olberg ||l99lj). This relati ve to the optical emission (Noriega-Crespo et al. 



jet presents a morphology in Ha and [SII] images that is |l996 
quite different from the morphologies of other well known 

stellar jets. While the best studied jets (e.g., HH 111 and Another special feature of the HH 110 jet is the lack of 
HH 34) show chains of aligned knots with more or less or- detection of any stellar source (suitable for powering the 
ganized, arc-like shapes, the HH 110 jet has a more chaotic J et ) alon S the outfiow axis - The morphology of HH 110 
structure: it starts in a well collimated chain of knots and first suggested that the driving jet source would be em- 
then widens in a cone of opening angle ~ 10°, with a rather bedded m a dark lane located to the north of the a P ex 
chaotic knot structure and appreciable wiggles along the of the outflow. However, searches at optical, near infrared 
length of the jet. The observed structure is reminiscent of and radio continuum wavelengths have failed to detect the 
a turbulent outflow. HH 110 drivin S source - 

Reipurth, Raga & Heathcote l|l996f) report the discov- 

Send offprint requests to: R. Lopez er y of another fainter jet, HH 270, located 3' northeast of 

* Partially based on observations made with the 2.5 m Isaac HH 110 - An embedded near-infrared source very close to 

Newton and the Nordic Optical Telescopes operated at the IRAS 05489+0256 was also found along the HH 270 flow 

Observatorio del Roque de los Muchachos of the Instituto de axis. They suggest that HH 110 is the result of the deflec- 

Astrofi'sica de Canarias. tion of the HH 270 jet (which travels in an E-W direction) 
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through a collision with a dense molecular cloud core. In 
this scenario, HH 270 strikes the molecular core and is 
deflected, giving rise to the HH 110 jet. The H 2 emission 
would then be a tracer of the region where the atomic jet 
an d the molecular cloud core interact (Noriega-Crespo et 
al. ll99ri) . 

Because of the interesting dynamics resulting from the 
interaction between a radiative jet and a dense obstacle, 
a number of theoretical papers on the H H 270 /HH 110 
system have been written. Raga & Canto l(l995h explored 
the initial stages of a jet/cloud collision with both analytic 
models and time-dependent, 2D numerical simulations. 
The final, steady state in which the jet has bored a hole 
through the (stratified) cloud w as stu died analytically and 
numer ically by Canto & Raga l)l996(l and Raga & Canto 
(1996). The initial state of a jet/cloud interaction (which 
we now appear to be seeing in the HH 270/HH 110 sys- 
tem) was re-explored, no w with 3D numerical simulations, 
by de G ouvei a dal Pino |l999). Hurka, Schmid-Burgk & 
Hardee (1999) carried out 3D numerical simulations of a 
jet interacting with a magnetized cloud. The most recent 
paper on models of a radiati ve jet /dense cloud interac- 
tion is the one of Rag a et al. £2002), who carried out 3D 
simulations with a simplified atomic/molecular network in 
order to obtain predictions of atomic and molecular emis- 
sion line maps for carrying out direct comparisons with 
the HH 270/HH 110 system. 

The HH 110 jet seems to be appropriate for search- 
ing the observational signatures of entrainment and tur- 
bulence through spat ially r e solved kinematic studies. 
Recently, Riera et al. ( 20oH l2003b|) presented detailed 
radial velocity studies obtained with multi-long-slit spec- 
troscopy and Fabry-Perot interferometry, and compared 
the observations with simple, para metriz ed models of a 
turbulent jet ( Canto, Raga & Riera 12003]) . 

In order to complement these radial velocity stud- 
ies, we now present new proper motion measurements of 
HH 110. These proper motions have been determined from 
deep [SII] CCD images that were obtained with a time 
baseline of 5424 days (i.e., nearly fifteen years), and show 
the kinematical properties of a large number of knots. The 
only previous proper motion determinations of HH 110 are 
those of Reipurth et al. I|1996l) . who study images with a 
shorter time baseline (~ 6 yr) and are not able to reach 
a level of accuracy and detail comparable to that of our 
present measurements. 

The paper is organized as follows. The observations 
are described in Sect. 2. In Sect. 3 we present the proce- 
dures to measure the proper motion velocities (Sect. 3.2) 
and discuss the obtained results (Sect. 3.1). The radial 
and total velocities are discussed in Sect. 4. Finally, the 
conclusions are given in Sect. 5. 

2. Observations 

In order to determine proper motions of HH 110 jet, a 
set of five CCD images has been used, yielding a time 
baseline of 5424 days (nearly fifteen years). Details of 



the epochs, telescopes, spatial resolution and exposure 
times are given in Tabled All the images were obtained 
through [SII] narrow-band filters, which included the [SII] 
A = 6717,6731 A emission lines. It is not possible to 
perform similar proper motion determinations from the 
Ha emission, as we did not obtain images through an Ha 
narrow-band filter during the observing runs carried out 
in the 1993 and 2002 epochs. 

The 1987 image was obtained at the Danish 1.5 m 
Telescope at La Silla. The 1988 image was obtained at 
the ESO 3.6 m Telescope. The 1994 image was obtained 
with EMMI at the ESO 3.5 m New Technology Telescope 
(NTT). Details on the acquisition and treatment of all 
these d ata ca n be found in the papers of Reipurth & 
Olberg l|l99l|) and Reipurth et al. l|l996|) . 

The 1993 image was obtained at the prime focus of the 
2.5-m Isaac Newton Telescope (INT) of the Observatorio 
del Roque de los Muchachos (ORM, La Palma, Spain). 
The detector used was a CCD with a coated EEV chip. 
A filter of central wavelength A = 6730 A and bandpass 
AA = 48 A was used. 

The 2002 image was obtained with the 2.6-m Nordic 
Optical Telescope (NOT) of the ORM. The Andalucia 
Faint Object Spectrograph and Camera (ALFOSC) was 
used, with a Ford-Loral CCD, and a filter of central wave- 
length A = 6724 A and bandpass AA = 50 A . Both the 
1993 and the 2002 images were obtained by combining five 
frames of 1800 s exposure each. For these two epochs, the 
individual exposures were reduced (bias-subtracted and 
flat-fielded) using the same procedure, through the stan- 
dard tasks of the IRAF 1 reduction package. The five indi- 
vidual frames were recentered using the position of several 
field stars, in order to correct for misalignments among the 
individual frames. Then the frames were median-averaged 
using the IMCOMBINE task of IRAF (with appropriate 
options to remove cosmic ray events). We obtained one 
deep [SII] image with a total integration time of 2.5 h for 
each of these two epochs (i.e., 1993 and 2002). These final 
images are not flux-calibrated. 

The central radial velocities of the main knots of 
HH 110 w ere obt ained from Ha Fabry-Perot data (see 
Riera et al. l2QQ3bl for more details of the observations and 
the data reduction). 

3. Proper motion measurements 

3.1. Method and results 

Proper motions of the main HH 110 knots have been de- 
termined from the CCD images listed in Tabled First, 
the five images were converted onto a common reference 
system and rebinned to the same pixel scale. We used the 
position of fifteen common field stars, in order to regis- 
ter the images. These reference stars are well distributed 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Reseach in Astronomy, Inc., under cooperative 
with the National Science Foundation 
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Table 1. Log of the Data 



Epoch 


Telescope 


Pixel scale 


Exp. Time 


Reference 






(arcsec) 


00 




1987 Dec 18 


1.5 m La Silla 


0.35 


3600 


1 


1988 Mar 5 


3.6 m ESO 


0.35 


1800 


1 


1993 Dec 16 


2.5 m INT 


0.55 


9000 


2 


1994 Jan 15 


3.5 m NTT 


0.35 


1200 


1 


2002 Oct 24 


2.6 m NOT 


0.19 


9000 


3 



References: (1) Reipurth et al. dl996l) : (2) Riera et al. <2003af) : 
(3) obtained by G. Gomez through the NOT Service Time facility. 
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20 40 

X offset (arcsec) 

Fig. 1. [SII] 6717+6731 image of HH 110 from the October 
2002 NOT image. Arrows indicate the proper motion ve- 
locity of each knot. Ellipses at the end of each arrow in- 
dicate the uncertainty in the components of the velocity 
vector. The scale of the arrows is indicated by the double 
headed arrow at the top right corner of the map, and cor- 
responds to a velocity of 100 km s _1 . Boxes indicate the 
regions used for the proper motion calculations. Offsets 
are measured from the peak position of knot A. North is 
up and East is to the left. 



around the HH 110 jet. The GEOMAP and GEOTRAN 
tasks of IRAF were applied to perform a linear transfor- 
mation with six free parameters that take into account 
relative translation, rotation and magnifications between 



Table 2. HH 110 Proper Motions' 



xvnot 


LJnser 


V T 


PA 


APA ° 




(arcsec) 


(km s _1 ) 


(de 


5) 


(deg) 


A 


0.0 


181.9 ± 5.2 


215 ± 


2 


21 


B 


7.9 


24.2 ± 3.7 


134 ± 


9 


-58 


13' 


14.9 


25.6 ± 2.6 


273 ± 


6 


79 


c 


22.8 


163.0 ± 5.8 


222 ± 


2 


28 


E 


40.1 


75.4 ± 4.0 


213 ± 


3 


19 


tr 


42.8 


153.5 ± 7.6 


206 ± 


3 


12 


H 


48.4 


133.3 ± 3.9 


220 ± 


2 


26 


Ie 


61.0 


44.4 ± 2.0 


197 ± 


1 


3 


I\v 


62.3 


144.7 ±22.3 


230 ± 


7 


36 


■In 


62.4 


159.8 ±17.8 


186 ± 


4 


-8 


Js 


70.0 


157.0 ± 7.3 


202 ± 


2 


8 


L 


73.0 


119.0 ± 7.7 


206 ± 


6 


12 


M 


81.3 


106.6 ± 2.9 


207 ± 


3 


13 


N 


91.9 


148.1 ± 4.0 


201 ± 


4 


7 


O 


101.8 


183.7 ± 5.4 


209 ± 


2 


15 


P 


113.6 


91.5 ± 2.7 


184 ± 


4 


-10 


Q 


121.2 


104.1 ± 5.2 


185 ± 


1 


-9 


R 


152.0 


84.9 ± 24.8 


128 ± 20 


-66 


S 


179.1 


68.9 ± 2.5 


216 ± 


2 


22 



a A distance of 460 pc has been adopted 

Offset from knot A (at the 1987 epoch) along the jet outflow 
axis 

c APA is defined as PA pm -PAj 



the frames. The final, transformed frames have a pixel size 
of 0"35. 

In order to check the accuracy obtained for the com- 
mon reference system of the five final images, two tests 
were performed: 

— Firstly, the displacements of the reference field stars 
in all images were calculated by computing the peak 
of the two-dimensional cross-correlation function be- 
tween pairs of frames (i.e., the first-epoch, 1987, and 
each of the other images) over small boxes around 
the reference stars. The average and rms of the dis- 
placements for the four pairs of images were less than 
0.09 ± 0.18 pixels for the x coordinate, and less than 
0.04 ± 0.14 pixels for the y coordinate. 

— Secondly, to test the accuracy of the transformations 
described above, we carried out astrometry for each 
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Fig. 2. Proper motions of the HH 110 knots. The offsets in the x direction are shown in the left panels and those in 
the y direction in the right panels. The knot is indicated in the top right corner of each right panel. For each epoch, 
a circle represents the displacement (in arcsec) in the ^-direction (left) and y-direction (right) of the corresponding 
knot, measured from the first-epoch image (which defines the origin of the time scale, set on 1987 Dec 18). Errors 
are indicated by the vertical bars. The least squares fits derived for each knot are shown by the continuous lines. The 
proper motions of the knots in the x, y directions are determined from the slopes of these lines. 

of the final images, using the (a, 6) coordinates of Catalogue 2 . With this procedure we confirmed that 
five bright field stars obtained from the USNO-B1.0 the pixel size indeed corresponds to 0"35 (see above) 



2 The USNOFS Image and Catalogue Archive is operated by 
the United States Naval Observatory, Flagstaff Station. 
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and that the orientation of each of the final images 
differs from the true North direction by less than 0?4. 

Thus, we concluded that the five final images used for 
the proper motion determinations were properly scaled 
and aligned, and were therefore suitable for estimating 
the proper motions of the HH 110 knots. 

In order to calculate proper motions for the knots of 
HH 110, we first defined boxes in each frame that included 
the emission from the individual knots (see Table [5] and 
Fig. ^ f° r the nomenclature of the knots and the defi- 
nition of the boxes). Then, for the four pairs of frames 
(i.e., the 1987 first epoch and each of the other images) 
we computed the two-dimensional cross-correlation func- 
tion of the emission within the pairs of boxes, and we 
determined the displacement in the x and y coordinates 
through a parabolic fit to the peak of the cross-correl ation 
function (see Reipurth et al. Il996l and Lopez et al. 1996 
for more details of this procedure). For each knot, the 
uncertainty in the position of the correlation peak was 
estimated through the scatter of the correlation peak po- 
sitions corresponding to boxes shifted from the nominal 
box by or ±2 pixels (i.e., 0"7) on each of its four sides 
(giving a total of 3 4 different boxes for each knot). We 
adopted as error for each coordinate twice the rms devia- 
tion of the positions of the cross-correlation peaks, added 
quadratically to the rms alignment error. 

The errors in the displacements obtained for the 1988 
image were much greater than for the rest of the images. 
This is probably a consequence of the poor quality of the 
image. Thus, we decided to discard the data obtained from 
this image from the rest of the analysis and use only the 
displacements from the 1987 image to the images of 1993, 
1994, and 2002. 

Finally, for each knot, we performed a linear regression 
fit to the displacements in each coordinate as a function of 
epoch offset from the first epoch (1987). For each coordi- 
nate of each knot we fitted a straight line to four points: a 
zero displacement for epoch offset zero, and the three mea- 
sured displacements. The fit was performed taking into 
account the error associated with each displacement. We 
assigned to the zero displacement an error equal to the 
quadratic average error of the three measured displace- 
ments. The fits obtained for the (x,y) coordinates of the 
different knots are shown in Fig. [21 From the slope of the 
fitted straight line and its uncertainty we determined the 
proper motion velocity (Vt) and position angle (PA) of 
each knot and their uncertainties. These values are listed 
in Table and displayed in Fig.Q A distance to HH 110 
of 460 pc has been adopted. 



3.2. Discussion 

Let us first define the HH 110 outflow axis passing through 
knots A to C, which define a well aligned structure. The 
outflow axis defined in this way has PAj = 193P6. From 
Figure ^ h is clear that, as a general trend, the proper 
motion velocities show a small westward deviation rela- 
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Fig. 3. Deviation of the proper motion velocity direction 
(PA pm ) from the axis of the HH 110 jet (PAj = 193?6, ref- 
erence dashed line ) as a function of distance to knot A. 
The dotted line gives the direction of the axis of HH 270 
relative to the HH 110 jet axis. Knots B and B' are ba- 
sically stationary (see text) and have been removed from 
the Figure. 



tive to this axis. This westward deviation tends to dimin- 
ish as one goes southwards from knot A. If on e com pares 
our Figure with Figure 4 of Reipurth et al. fl996), one 
sees a general agreement between the two sets of results. 
However, there are substantial differences for some of the 
knots, which are discussed in detail below. 

At distances > 130" south from knot A the emitting 
knots show large deviations with respect to the outflow 
axis, with a locus that first curves to the E and then to 
the W, becoming again more or less parallel to the out- 
flow axis. In this region, we have measured the proper 
motions of knot R (which is moving to the SE, in a di- 
rection parallel to the locus of the jet) and condensation 
S (again moving parallel to the locus of the jet, which in 
this position is also approximately parallel to the outflow 
axis). The resulting structure is quite striking, as it ap- 
parently shows that the knots are flowing along a curved 
"channel" , in contrast to the expected behaviour of a pure 
ballistic motion (note that the knot proper motions are 
always parallel to the local direction of the locus of the 
emitting region). However, as can be seen from Figure ^ 
this conclusion is mostly based on the SE proper motion 
of knot R, which has a highly uncertain measured proper 
motion. If we remove this knot from the analysis, we would 
only conclude that all of the knots of HH 110 have proper 
motions directed sligthly to the W of the outflow axis. 

In order to illustrate the westward shift, we calculate 
the differences APA=PA pm — PAj between the directions 
of the proper motions and the direction of the HH 110 axis, 
PAj = 193°. 6 (see above). In Figure[j|]we plot these angu- 
lar differences as a function of distance along the outflow 
axis. In this figure we see that knots P, Q and Jn move in 
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Fig. 4. Moduli of the proper motion velocities of the 
HH 110 knots, as a function of distance to knot A. 



a direction to the E of the HH 110 axis. Knot R also moves 
to the E (but, as we have discussed above, it is not clear 
whether or not this proper motion measurement is signif- 
icant), and the remaining 13 knots with measured proper 
motions move in a direction to the W of the HH 110 axis. 
In this figure the corresponding directions of the HH 110 
and HH 270 jet axes have been indicated by dashed (for 
HH 110) and dotted (for HH 270) lines, respectively. The 
directions of the motions of these 13 knots lie in between 
the HH 110 axis and the HH 270 axis (the jet that is 
apparently deflected through a collision in order to form 
HH 110). In addition, this plot shows that while for the 
knots with distances d < 50" from knot A the proper mo- 
tions have PA pm — PAj ~ 10° to 30° (with an average value 
(PA pm - PAj) = 21°), for distances 50" < d < 130" the 



-10° to 36° (with an 
This decrease in the 



proper motions have PA pm — PAj ~ 
average value (PA pm — PAj) = 7°) 
average angular offset (PA pm — PAj) between the proper 
motions and the outflow axis shows that there is a partial 
convergence between the proper motions and the outflow 
axis as one moves along the HH 110 jet. However, at dis- 
tances d > 130" from knot A we have knots R and S (see 
Figure |3J) , that again have proper motions that are mis- 
aligned with the outflow axis. These proper motions, as 
stated above, appear to be aligned with the local direc- 
tion of the curved locus of the jet. 

In ord er to model the HH 270/110 system, Raga et 
al. I|2002(l performed 3-D gasdynamics simulations for a 
radiative jet/cloud collision. Two assumptions for the in- 
cident jet velocity were considered: a time- independent, 
constant direction velocity (Model A) and a precessing jet 
of sinusoidally varying velocity (Model B). In both mod- 
els, the proper motion velocity vectors deviate from the jet 
axis direction, and the predicted deviation trend success- 
fully reproduces the observed behaviour in HH 110. Thus, 
it might be reasonable to interpret the observed lack of 
symmetry by a jet/cloud collision. 



In Figure 21 we show the magnitude of the proper mo- 
tion velocities as a function of distance d from knot A. In 
the d < 90" region, there are 4 low velocity knots (knots 
B, B', E and Ie)- All of the other knots appear to fall on 
a general trend of decreasing velocities as a function of d 
(with knot A having Vt — 180 km s _1 and knot M, at 
d = 81"3, having V T ~ 110 km s" 1 ). For d > 90", the 
proper motions of the knots have a wide range of values, 
from ~ 70 km s _1 for knot S up to ~ 18 5 km s -1 for knot 
O. Note that for knot O Reipurth et al. l|l99fil) determined 
a velocity of only 105 km s . 

Knot B seems have remained nearly stationary 
through the fifteen years elapsed between the first-epoch 
and the last-epoch frames, at least within the errors of 
our proper motion determinations, as shown in Figure Q 
and Table El In order to explore the "stationarity" of this 
region of the jet, we marked as B' the emission located 
between knot B and knot C, and we evaluated a proper 
motion for B'. The proper motion velocity found for B' 
is similar to the velocity of knot B. Thus, it appears that 
HH 110 has a region of a length of ~ 14" along the jet 
direction ( i.e., the region going from the end of knot A to 
the beginning of knot C) without appreciable proper mo- 
tions between the first-epoch and last-epoch optical im- 
ages. The lack of motion of k not B is in agreement with 
the results of Reipurth et al. l|l99fil) . 

A possible explanation for the nature of the region 
around knot B is that it might correspond to the location 
of the jet collision with the cloud. In fact, it should be 
noted that the values of the [Nil] /Ha emission line ratio 
derived from long-slit spectra around knot B are higher 
than the values found in its surrounding s, arou nd knots A 
and C (Reipurth et al. ll996t Rie ra et al. l2003al Figure 4). 
In addition, Riera et al. ( 2003a) found a minimum in the 
[SII]/Ha line ratio around knot B relative to its surround- 
ings. This trend in the [Nil] /Ha and [SIIJ/Ha line ratios is 
indicative of a higher gas excitation around knot B relative 
to its surroundings and gives support to this interpreta- 
tion, although alternative explanations are also possible 
(see next section). 

Knot I is elongated in the E-W direction. It has two 
peaks, which we label Ie and lw ■ The ratio of the relative 
intensities of these two maxima reverses from the 1993 to 
the 2002 images (with I^; being brighter in 1993 and lw 
being brighter in 2002). We find that while lw shows a 
proper motion similar to the ones of other nearby knots, 
Ie shows a much lower proper motion (see Figure 0J. 
Therefore the proper motion of I^; appears to be some- 
what doubtful, and we suspect that it might be the result 
of a strong change of morphology of the emission in this 
region, rather than a real motion of a physically well de- 
fined knot. 

Finally, for knot A we confirm the change in morphol- 
ogy with time pointed out by Reipurth et al. lll996h . In 
Figure [5] we show a superposition of the 1987, 1993 and 
2002 images of knot A. From this figure we see that the 
morphology of this region changes quite substantially with 
time. Not only does the peak of the emission move, but 
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Fig. 5. Close-up of HH 110 knot A from images of Dec. 
1987 (grey overlaid), Dec. 1993 (dashed line contours) and 
Oct. 2002 (solid line contours). North is up and East is to 
the left. 

also knot A develops a "tongue" of emitting material ex- 
tending along the outflow axis. We obtain a proper mo- 
tion of ~ 180 km s _1 for knot A, which is highe r tha n 
the ~ 50 km s" 1 determined by Reipurth et al. ([[996). 
Despite this change in the morphology of knot A, from 
Figure it is evident that the peak of knot A shifts as a 
function of time. This displacement of the peak produces 
the high proper motion which we find for knot A. 

4. Radial and full spatial velocities 

From the long-slit spectra obtained along the HH 110 jet 
axis and across several of its brightest knots, we found 
that the values of the radial velocity derived from the Ha 
emission differ by less than 15-20 km s _1 (for a spec- 
tral resolution of ~ 20 km s^ 1 ) from the values derived 
from the [SII] emission (see, eg, Riera et al. l2003al Figure 
3). Thus, it seems reasonable to combine proper motions 
derived from [SII] imaging and radial velocities derived 
from the Ha emission line to evaluate full spatial veloci- 
ties of the knots. Accordingly, we used the more complete 
HH 110 Fab ry- Pero t data of the Ha emission line from 
Riera et al. ( 2003b) to obtain the heliocentric radial ve- 
locity of the HH 110 knots. ^From these data we computed 
the line profile integrated over the boxes defined for mea- 
suring the proper motions of the knots (see Figure 0, and 
determined the radial velocity of the emission peak with 
a parabolic fit. The results are shown in Figure where 
we present the radial velocities with respect to the sur- 
rounding cloud (which has a heliocentric radial velocity of 
+23 km s _1 ) as a function of distance from knot A. 

We can clearly identify two groups of knots in Figure|Hl 
The first group includes knots A, B, C and E, which show 
an increase of the value of the radial velocity with dis- 
tance (from ~ —37 km s _1 at knot A to ~ —24 km s _1 
at knot E). The second group corresponds to the knots 
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Fig. 6. Radial velocities of HH 110 knots, as a function of 
distance to knot A. The typical value of the velocity error 
is ~ 10 km s _1 . 

found at distances from 50" to 130" from knot A (includ- 
ing the knots located between knot H and knot Q), which 
show radial velocities ranging from —38 to —53 km s _1 . 
There is a trend of more negative radial velocities at larger 
distan ces, which has been previously reported by Riera 
et al. l|2003al:l2003b|) . Finally, knot S shows a velocity of 
~ —25 km s _1 relative to the surrounding cloud. 

From the radial velocities (relative to the surround- 
ing cloud, see above and Figure^) and the proper motion 
velocities (see Table FJJ we calculated the total spatial ve- 
locity v to t of the knots, and the angle <p between the knot 
motion and the plane of the sky (with positive values of <p 
towards the observer). The results of this estimation are 
shown in Figure From this figure we see that the ori- 
entation angle has a value (j> — 9° to 12° for the region 
within 30" from knot A. At larger values of d (> 40"), the 
orientation angle has values <fi — 15° to 25°. This range of 
~ 10° in orientation angle is consistent with an opening 
angle of ~ 10° (in t he plan e of the sky) of this region of 
HH 110 (Riera et al. l2003bh . 

A general trend of decreasing spatial velocity (i>tot) 
with distance (d) can be appreciated from Figure [7| The 
velocities decrease from w to t — 190 km s^ 1 for knot A to 
t>tot — 70 km s _1 for knot S. Knots E and O are the ones 
that show the largest deviations from the general trend of 
decreasing velocities as a function of distance from knot 
A. 

Regarding knot B, note that the region around this 
knot is not fully stationary considering the 3-D spatial 
velocity, since we determined a continuous increase of the 
radial velocity from knot A to knot C. Thus, instead of 
being the location of the jet collision with the cloud, an 
alternative speculative explanation might be that the jet 
is excavating a channel around knot B, with its entrance 
around knot A, coming out around knot C. 

It should be noted that we carried out the evaluation 
of the spatial velocity vtot of the knots by assuming that 
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Fig. 7. Total velocities (upper panel) and orientation an- 
gle, denned by the direction of the knot motion and the 
plane of the sky (lower panel), as a function of distance 
to knot A. The typical value of the total velocity error is 
~ 12 km s _1 . The typical value of the orientation angle 
error is ~ 6°. 

v p , the pattern velocity of the knots (i.e., the knot proper 
motion along the jet axis, v v = Vr/ cos <p) is similar to the 
fluid velocity, Vf. This condition is fulfilled in the case of 
highly supersonic, dense jets, where the knots represent 
internal working surfaces (IWS). However, if the knots 
originate from Kelvin-Helmholtz (K-H) instabilities, then 
the ratio rj — Vf/v p could be higher than unity (eg, up 
to a factor 2, Micono et al. Il998|) . If the HH 110 knots 
originated from K-H instabilities instead of corresponding 
to IWS, the results we derived give a lower limit to the full 
spatial velocity i>tot> an< ^ the t rue jet motion direction </> 
would lie closer to the plane of the sky (i.e, given the low 
radial velocities found for the HH 110 knots, the values 
derived for v to t would have to be scaled by a factor of 
the order of 77, while the derived <p values would have to 
be scaled by a factor of the order of f]^ 1 ). This rj scaling 
value cannot be measured from the present observations, 
and has to be assumed from a model. 

The result we found for the behaviour of the spatial ve- 
locity («tot) is very int eresting , because it reverses the con- 
clusion of Riera et al. (2003afl. who deduced an increasing 
velocity as a function of distance along the HH 110 flow. 



The ori gin of this difference lies in the fact that Riera et al. 
(2003a) took the measured radial velocities and converted 
them into a full spatial velocity by deprojecting these ra- 
dial velocities assuming the same value of </> = 35° for the 
orientation angle (with respect to the plane of the sky) for 
all of the knots. The fact that we now deduce a position- 
dependent orientation angle for the knots along HH 110 
(see Figure 01 modifies the result of Riera et al. l|2003a|) in 
such a way that the observed velocity vs. position trend 
is now reversed. 

The observed trend in the behaviour of the total ve- 
locity is not predicted by the jet/cloud models of Raga 
et al. l)2002() . In particular, Model A gives velocity values 
significantly lower than the values measured in HH 110 
and, in addition, predicts an increase of the proper mo- 
tion velocity with distance from the cloud (i.e., with the 
distance from knot A). This prediction is the opposite of 
the observed trend (i.e., decreasing proper motion veloci- 
ties away from knot A). 



5. Conclusions 

We carried out proper motion measurements of the knots 
of HH 110 on a series of five [SII] 6717+6731 CCD images 
covering a time baseline of ~ 15 years. With these mea- 
surements we obtained a better accuracy and were able 
to describe the motion of more knots than in the p reviou s 
proper motion determinations of Reipurth et al. (1996). 
^From the Fabry-Perot data of Riera et al. I|2003b|) . we 
obtained the average radial velocities for the boxes around 
the knots that were used for determining the proper mo- 
tions. Combining these radial velocities with the proper 
motions, we obtained the full spatial velocities «tot and 
the orientation angles (f> (with respect to the plane of the 
sky) of the different clumps. 

We found a number of interesting effects: 

— The proper motions of the knots in the northern part of 
the outflow (i. e., close to knot A) have an orientation 
angle in between the HH 110 and the HH 270 axes, 
westward of the HH 110 axis. 

— The knots further away from knot A have proper mo- 
tions which are better aligned with the HH 110 axis. 

— The proper motions show a strong lack of symmetry 
with respect to the outflow axis, and show a general 
tendency to point towards the west. This trend is pre- 
dicted by the jet/cloud collision models of Raga et al. 
l|2002h . 

— There is a region of ~ 14" in length around knot B 
without appreciable proper motion velocity. We spec- 
ulate that this region might trace the location of the 
jet collision with the cloud. Signatures of a higher gas 
excitation around knot B relative to its surroundings 
are found from the [Nil] /Ha and [SII] /Ha line ratios, 
which is in favour of this hypothesis, although other 
explanations are possible. 

— Knot R, which is located in a region in which the jet 
has a highly curved shape, appears to be aligned with 
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the locus of the jet (deviating substantially from the 
direction of the outflow axis). This result is somewhat 
marginal, because the proper motion determined for 
knot R is quite uncertain. 
— The spatial velocity v to t shows a general trend of de- 
creasing velocities as a function of distance from knot 
A (i.e., accounting correctly for the changes in orienta- 
tion angle, one finds a "deceleration" of the full spatial 
velocity along the HH 110 flow). This trend in the be- 
haviour of the total velocity is not pre dicted by the 
jet/cloud collision models of Raga et al. l)2002|) . 

We conclude that HH 110 is a jet with somewhat pe- 
culiar characteristics. It should be noted that the spatial 
velocities of the knots along HH 110 show a decrease by a 
factor of ~ 3 between knot A and knot S. This slowdown 
occurs over a distance of ~ 10 18 cm, which is much shorter 
than the distances over which one expects the braking 
due to interaction with the surro unding environment to 
becom e important (Cabrit & Raga l2002t Masciadri et al. 
12002ft . Therefore, the strong slowdown of the spatial veloc- 
ity along HH 110 might be an evidence for an anomalously 
strong interaction between the outflow and the surround- 
ing environment (such as might occur, e. g., during the 
deflection in a jet/cloud collision). Alternatively, the ob- 
served slowdown could be the result of variability in the 
ejection velocity, or of changes in the jet/cloud impact 
region. In order to discern between these different possi- 
bilities, it will be necessary to carry out further modelling 
of the HH 270/110 flow. 
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